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Abstract 



CN ' We will numerically investigate two-particle correlation function of CERN- 



SPS 158 A GeV Pb+Pb central collisions in detail based on a (3+1)- 
CH ■ dimensional relativistic hydrodynamical model with first order phase tran- 

^rr \ sition. We use the Yano-Koonin-Podgoretskii parametrization as well as the 

<^ \ usual Cartesian parametrization and analyze the pair momentum dependence 

of HBT radii extracted from the parametrizations. We find that the inter- 
pretation of the temporal radius parameters as the time duration in YKP 
parametrization is not available for the hydrodynamical model where the 
source became opaque naturally because of expansion and surface dominant 
freeze-out. Finally, effect of the phase transition on the source opacity is also 
discussed. 
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I. INTRODUCTION 

Ultra-relativistic heavy ion collision aimed at producing the Quark-Gluon Plasma (QGP) 
state is one of the most attracting problems in modern nuclear-particle physics 0. In these 
decades, many kinds of candidates for the signal of the production or existence of QGP are 
proposed. Two-particle correlation has been accepted as one of the probable observables for 
finding anomalous events in which QGP can be produced, because of the naive speculation 
that an extremely huge hot region will be produced if phase transition occurs @Q|. In 
quantum optics, estimating source size through quantum correlation of emitted particles is 
a well-established method called the Hanbury-Brown Twiss(HBT) effect H|J. 

Usually, parameters obtained from the Gaussian fit of the correlation functions are re- 
garded as a radius of the particle source (HBT radius). However, if the source is not static, 
the interpretation of the Gaussian parameter is not straightforward @J^]. In the case of 
ultra-relativistic heavy ion collision, in addition to the highly dynamical properties, many 
kinds of experimental limitations, such as rapidity windows, make the situation more com- 
plicated. In the recent analyses, there exist two ways of parametrization in the Gaussian fit 
of the two-particle correlation function; one is Cartesian parametrization MM , and the other 



is Yano-Koonin-Podgoretskh parametrization |10|,|Tl],|T2[ . In both parametrizations, we can 
obtain their respective "size" of particle sources. However different parametrization gives 
different values and interpretation of their values is not clear enough. 

In this paper we evaluate two-particle correlation functions of pion based on a hydro- 
dynamical model with phase transition. Parameters in the hydrodynamical model are so 
tuned as to reproduce single particle distributions of CERN 158 A GeV Pb+Pb collisions 
T3|j . Comparing the source size obtained through correlation and freeze-out hypersurface 



of the hydrodynamical model, we investigate the physical meaning of the HBT radii in 
both parametrizations. The opaque property ]14[ of the source is also analyzed based on 
a hydrodynamical model. In a hydrodynamical model, particle emission takes place on the 
freeze-out hypersurface which is thin but not an instant in time. Expansion of the fluid also 
should naturally occur. Hence, the apparent opaque property of the freeze-out hypersurface 
is not trivial. The influence of first order phase transition is also investigated. 

In section |T[ we make a brief introduction to the hydrodynamical model with phase 
transition. In section IDT], we formulate two-particle correlation functions. Transverse mo- 
mentum dependence and rapidity dependence of source parameters are discussed in Section 
T\/\ and |V|, respectively. Section [Vl| is devoted to the concluding remarks. 



II. HYDRODYNAMICAL MODEL 

Focusing our discussion to the central collisions, we may assume cylindrical symmetry to 
the system and introduce cylindrical coordinates r, rj, r and <ft instead of the usual Cartesian 
coordinates x, y, z and t, where, 

t = r cosh rj, (la) 

z = rsinh?7, (lb) 

x = r cos0, (lc) 



y = r sin0. (Id) 

The four velocities of the fluid U^ is given as, 

U* = U T cosh 7] + U v sinh rj, (2a) 

U z = U T sinh n + U v cosh 77, (2b) 

U x = U r coa<l> + U*sm<l), (2c) 

£^ = £/ r sin0-£^cos0. (2d) 

Taking the constraint, U^U^ = 1, and cylindrical symmetry into account, U^ can be 
represented with two variables YL and Ft |T5| , 



U T = cosh Yr cosh(F L — ??), (3a) 

If = cosh F T sinh(F L - 77) , (3b) 

U r = smhY Tl (3c) 

C/* = 0. (3d) 

As is well known, the relativistic hydrodynamical equation is given as 

d,T^ = 0, (4) 

and the energy momentum tensor of perfect fluid is given by 

T^ = EU»U U - P{g^ - U»U U ) (5) 

where E and P are energy density and pressure, respectively, and the metric tensor g^ u is 
g/iv _ diag.(l, — 1, — 1, — 1). In a hydrodynamical model, thermodynamical quantities such 
as E and P are treated as local quantities through temperature T = T(x M ). 

In order to solve the above hydrodynamical equation, we need an equation of state. 
In this paper, we adopt the bag model equation of state as the simplest model for the first 
order QCD phase transition. In the bag model, thermodynamical quantities below the phase 
transition temperature correspond to a massive hadronic gas. On the other hand, in the 
high temperature phase, thermodynamical quantities are given by a massless QGP gas with 
bag constant B. In this paper, we assume pions and kaons in the hadronic phase and u-, d-, 
s-quarks and gluons in the QGP phase. Putting the phase transition temperature, T c = 160 
MeV, into the condition of pressure continuity, we can obtain the bag constant as B = 412 
MeV/fm 3 . 

At the phase transition temperature T c , pressure P is continuous but other quantities 
such as energy density E and entropy density S have discontinuity as a function of temper- 
ature. We parameterize these quantities during the phase transition by using the volume 
fraction A [|I(|,|17]]. We assume that thermodynamical quantities at T = T c are given as the 



functions of the space-time point through \{x 



A' 



E(X) = XE QGP (T C ) + (1 - A)£ H ad(T c ), (6a) 

S(X) = XS QGP (T C ) + (1 - X)S HAD (T C ), (6b) 



where < A(x M ) < 1. This parameterization enables us to solve the hydrodynamical equation 
easily. 

We consider that the hadronization process occurs at freeze-out temperature T{ = 140 
MeV. In the evaluation of the one-particle distribution in a hydrodynamical model, the 



Cooper-Frye formula [[HI] is often used, however, there exist several delicate matters on 



the treatment of freeze-out hypersurface based on this formula [|19|,^(J. In this paper, we 
analyze not only one-particle distribution but also two-particle distribution, hence, we adopt 
another type of the simplest formula which enables us to extend two-particle distributions 
straightforward, 

d 3 ANi _ U,jf 1 U x dcr x 



d P 3 ^p 2 + m 4 2 exp(^r)-l (2 



7T) 3 



where p° = Jp 2 + mf and the subscript i stands for pion or kaon. Integrating (0) on the 
freeze-out hypersurface E, we obtain the rapidity distribution, 

dNi r U^ U x da x p T dp T d(j) 



dY h eX p (Zf.) - 1 (2 



m 



,3 



and the transverse mass distribution, 

1 dNi f Utfp U x da x dYd<p 



m T dm T 7e exp ( ^f- ] - 1 (2tt) 



T 



(9) 



In a previous paper, we have already analyzed one-particle distribution of hadrons by using 
results of our hydrodynamical model and compared it to the data at CERN SPS experiment 
]T3| . The detailed structure of space-time evolution and particle distribution is discussed in 
the Appendix. 

III. TWO-PION CORRELATION FUNCTION 

Two-particle correlation function is defined as 

c(k k) ^(kx,k a ) 

C ^*>) - w{kl)w{k2 y ( 10 ) 

where W(k) and W(ki, k 2 ) are one-particle distribution E k dN/d 3 k and two-particle distri- 
bution E kl E k2 dN/d 3 k 1 d 3 k 2 , respectively. Assuming that the source is completely chaotic, 
for simplicity, the two-particle correlation function for identical bosons (i.e.,pions) is rewrit- 
ten as plj 

where 



J(k!,k 2 ) = f U fl da^e^ kl - k2) ^^/F{k 1 ,x)F{k 2 ,x) 



(12) 



and 



F(k,x) 



U^ 



1 



(2vr) 3 E k exp 



Uvky_ 



(13) 



T 



f 



As usual, we introduce the average momentum K^ = \{k^ + k^) and the relative mo- 
mentum g M = ki — k% of two particles. Taking an approximation that ki ~ k 2 — K in 



Eq.flllD, the two-particle correlation function C(k 1; k 2 ) can be expressed as J8|j22 

C(g", #") ~ 1 + 

where S(x^, K^) is the source function 



Jd 4 xe iq ^ S{x^,K^ 
\fd*xS{x».K»)\ 2 



S(x^K^) = f U tl {x')da tl {x')E K F{K,x')8 i {x 



x 



(14) 



(15) 



Furthermore, expanding e l<?f ' xM in Eq.(|l4]) with q^x^ <C 1, we can obtain the simple expression 

C{<f,K») - 1 « 1 + (q^) 2 - ((q^) 2 ), (16) 

with (A(x (i )) being weighted average , 

Sd i xA(x tt )S{x' l ,K' 1 ) 



(A{x*)) 



Jd 4 xS(x»,K») 



(17) 



Consequently, we can recognize the two-particle correlation function as the second order mo- 
ment of g M x M with the weight of S(x ,x , K^). Because of the on-shell property of final particles, 
the relative momentum q^ has three independent components. Choosing the three indepen- 
dent components of g M in an appropriate manner, we can obtain meaningful information on 
the space-time structure of the particle source. 

In the estimation of the source sizes from the two-particle correlation function (11), two 
types of Gaussian parametrizations \TT\ are widely used. One is the Cartesian parametriza- 
tion, including the outward-longitudinal cross term HI, and the other is so called Yano- 



Koonin-Podgoretskh parametrization |I0|,I2"|. We take the z-axis in the direction of the 
collision and the x-axis parallel to the transverse component of the average momentum 
K^ = (K°,Kt,0,Kl), as usual. Then the z-axis, x-axis and y-axis correspond to so-called 
"longitudinal", "outward" and "sideward" directions, respectively. 

In the Cartesian parametrization, the three independent components of the relative mo- 
mentum q^ are (q x ,q y ,qz) = (<Zout, fele, QWg)- In this parametrization, the fitting function 
with the average momentum K^ is given as 



C(q^K^) = 1 + A exp f-i4ie(^k 2 ide - *L(tf")?L - KnAK^qL, ~ 2R 2 ol (K^q ont q h 



<>ns 



(18) 



There are four fitting parameters in Eg. (|T8D , namely, -R s ;de, -Rout, -Riong and R Q \, which are 
so-called HBT radii. As for A, we fix its value one because we have assumed chaotic sources. 
In this parametrization, the meaning of these HBT radii, which is given through Eq. (|T6|) . 
depends on the observer's longitudinal frame. The coordinate system where $ = K^/K = 
is called longitudinal co-moving system (LCMS), and on this system, the interpretation 
becomes very simple l|, 



E>2 
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(19a) 
(19b) 
(19c) 
(19d) 



where x = x — (x) and /3± = Kt/K°. It is clearly seen that these radii mix spatial and tem- 
poral components. Therefore, from this parametrization it is not easy to extract accurately 



physical information such as emission duration |L2 



In the Yano-Koonin-Podgoretskh (YKP) parametrization, the three independent com- 

and q° 



ponents of relative momentum are chosen as q± 



The fitting function is given as |12 



Iside + Qout, Ql 



Qh 



ong 



E 1 — E 2 . 



C(q», K") = 1 + Aexp {-R 2 ± {K»)q 2 ± - R\(K») 



0\2 



it - W) 



R 2 (K») + R\(K»)][q^(K»)} 2 } 

(20) 



where u(K») = j(K^)(l, 0, 0, v(K»)), and j(K^) = l/yjl - v 2 (K^). We fixed A = 1 as well 
as in the Cartesian parametrization. Hence, the fitting parameters are R±, Ro, R\\ and the 
so-called Yano-Koonin-Podgoretskh velocity v which appears in the four-velocity w M . The 
first three parameters are the HBT radii that are invariant under longitudinal boosts. This 
is one of the most important properties of the YKP parametrization. The four- velocity u^ 
has a longitudinal component only, and the Yano-Koonin-Podgoretskh velocity v can be 
interpreted as the longitudinal source velocity JT5 1 . The interpretation of the HBT radii in 
the YKP parametrization can be given in the same manner as in Cartesian parametrization. 
Denoting deviation as (AA) 2 = ((A — (A)) 2 ) and 



A = (At) 2 - -f(xt) + -I [(Ax) 2 - (Ay) s 
Px Px 

B = (Az) 2 - ^l(xz) + % [(Ax) 2 - (Ay) 







0i 



0\\ 



C = (zt) - %L(tx) - -z-(zx) + ^ (Ax) 2 - (Ay) 











0i 



(21a) 
(21b) 
(21c) 



the HBT radii and Yano-Koonin-Podgoretskh velocity are expressed as follows [12 



v(K») 



A + B 
2C 



2C 
A + B 



(22a) 



rKk") 

R\{K^) 



(Ay) 2 

A-vC 

B-vC. 



(22b) 
(22c) 
(22d) 



Even in the YKP parametrization, R± directly corresponds to transverse size. In the Yano- 
Koonin-Podgoretskh(YKP) frame which is defined by v = 0, these expressions ( 22b - 22d" ) 
become very simple, 

(23a) 
(23b) 



(23c) 



The values of these quantities do not change because they are invariant under longitudinal 
boosts. But now the (• • •) is taken in the YKP frame. It has been shown that the last 
two terms in Eqs.( 23b" ) and ( |23cj ) can be considered as a small correction within some class 
of the thermal model |T2J. Hence, the three HBT radii can be considered to measure the 
"source size" directly. Specifically, Rq measures the emission duration on which, according 
to a naive expectation, the effect of first order phase transition may be observed. Of course, 
it is not clear that this interpretation works well even in the hydrodynamical model. 



Rl(K^) = 


= (Ay) 2 • 






Rl(K») - 


= (At) 2 - 


-it m+ wA iAx)2 - 


(Ay) 2 ; 


Rl(K") = 


= (A.) 2 - 


"^(-) + ^[(Ax) 2 


- (Ay) 2 ; 



IV. TRANSVERSE MOMENTUM DEPENDENCE OF SOURCE PARAMETERS 

In this section, we will discuss the dependence of HBT radii on the transverse momentum 
Krp. Because of the lack of experimental data with fixed average momenta, in order to make 
a comparison with the experimental data, we integrate the distribution functions in ( [TT]) 
with respect to average momenta and 



C(g"; K T ) 

C(q»;Y) 



1 + 



1 + 



/c^|I(k 1; k 2 )| 2 
JW L Vr(ki)W(k 2 )' 

S K T dK T W(\L X )W{\L 2 ) 



(24a) 
(24b) 



is used. 

Figure [I] shows the dependence of Cartesian HBT radii on the average transverse mo- 
mentum Kt- Closed circles stand for our numerical results obtained from the correlation 
function ( |24a| ) through the fitting function (|18D . Experimental results taken from NA49 
preliminary data |23) are plotted as open circles. Our results almost agree with the experi- 
mental data. However, we get the HBT radii a little larger than experimental results at high 
Kt- In the upper left figure of Fig. IB, though experimental -R s ;d e becomes smaller at larger 
Kt, this tendency appears only weakly in our simulation. This Kt dependence of -R S ide can 
be explained by the existence of strong transverse flow H and contributions from resonance 
decay |25j. Hence, the differences between our calculation and experimental results might 



be caused by the simplification that neglected the resonance decay. 



As we discussed in the previous section, a simple interpretation is not available to HBT 
radii in the Cartesian parametrization; investigating the source parameters based on the 
YKP parametrization is important. Figure shows the dependence of YKP radii on average 



transverse momentum Kt- As mentioned in Sec. |T|, in the YKP parametrization, the HBT 
radii are expected to stand for longitudinal invariant "source size" . In Fig. |2|, closed circles 
stand for YKP radii obtained from the correlation function ( p4a| ) through fitting based on 



Eq.(^Dj). Solid lines stand for the space-time extensions of freeze-out hypersurface evaluated 
from the right-hand side of Eq.( |23a|) -( |23c|) . These quantities are expected to agree with 
our HBT radii because in our calculation we consider only thermal pions jpfifl . Dotted lines 
shows "source sizes" (Az for R\\ and At for R Q ) calculated from the second order moment, 
i.e., first terms only in the right-hand side of Eq.( [23b|) - (|23c|) . In the calculation of space- 



time extensions and "source sizes" , average momenta are fixed at the central value of the 



integration range. NA49 experimental data p4j are plotted as open circles. 



The upper figure shows the dependence of the transverse source radii R± on Kt- As we 
have already discussed, R± in the YKP parametrization is the same as -R S id e m the Cartesian 
parametrization. The difference between the space-time extension and YKP radius R± is 
about 0.5fm. 

The middle part of Fig. [2] shows the dependence of longitudinal source radius R\\ on Kt- 
Our result from the fit is consistent with both the experiment and space-time extension. 
Furthermore, the "source size" Az also agrees with the space-time extension. Hence, the 
usual interpretation of R\\ in the YKP parametrization works well in our calculation. 

The lower figure shows the dependence of the temporal source parameter Rq on Kt- Our 
results from the fit are a few fm smaller than experimental results and become smaller as Kt 
increases. The space-time extension calculated from Eq. ( |23b|) shows similar behavior in Kt 
but quantitatively agreement is not so good at large K T . Because of the on-shell condition, 
the available g M region is limited and a great ambiguity in fitting process can occur. For 
example, R% sometimes tends to become negative if we take the most likely values without 
any restriction to Rq. However, it becomes small but positive definite values if we impose 
the restriction on R^ to be positive. As for the "time duration" At, though it agrees with 
R at low Kt, it becomes larger than R Q at large Kt- Hence, the interpretation that i?o 
measures the time duration directly in the YKP parametrization is not correct in this region. 



This can be explained by the effect of "source opacity" |T3J, which will be discussed later. 
On the other hand, the values of Rq and At are about 2-3 fm, obviously smaller than R±, 
Rn . Hence, we can say that clear evidence of the first order phase transition does not appear 
in the time duration. It was already pointed out that the large time duration (or longer 
lifetime) might not occur for Bjorken-type scaling, the initial condition we used [0,0- 

The "opaque" source emits particles dominantly from the surface and the "transparent" 
source emits from whole region. Whether the "real" source is opaque or transparent is an 



interesting problem P?jj] . Aiming for this point, we analyze the source function S(x^,K fl ) 
for Kt = 450 MeV and Y n7T = 4.15. Figure |3] shows the source function as a function of 
transverse coordinates x and y. Here the source function is integrated in the temporal and 
longitudinal coordinates, St{x,v) = * J dzdtS(x^,K^), and normalized on the transverse 
plane as / dxdy St{x,h) = 1. In this figure, emitted pions have only x component of 
momentum. Figure ^| clearly shows that pions are emitted mostly from the crescent surface 
region and the source is thinner in the x-direction than in the ^-direction. This is the typical 
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property of the opaque source. Therefore, the difference between Ax and Ay, that is, the 
difference of space-time extensions (Ax) 2 — (Ay) 2 is a good measure of the opaque property 
2JJ. 

Figure f| shows the dependence of the second and third terms in Eq . (|23b|) on K?. The 
third term ((Ax) 2 — (Ay) 2 )/ ' 0]_ has a large negative contribution to Rq at a large Kt region, 
while the contribution from the second term —2(xt)/j3± is not so large. Hence we may not 
consider the third term as a small correction to (At) 2 and this is the reason why the simple 
interpretation of Rq in the YKP parametrization is not available in our calculation. A simple 
model analysis ]2I| has shown the contradictory result that R diverges to — oo as K T — > 
(P± — » 0) for the opaque source. This is due to the tendency of the denominator of the third 
term in Eq. (|23b| ) to go to zero as Kt — * while the numerator is kept as finite. However, 
this divergence does not occur in our model because the numerator, (Ax) 2 — (Ay) 2 , goes to 
as K T -»• 0. 

Figure [| also shows the source function St(x, y) with the artificially neglected transverse 
flow. In order to see the effect of transverse flow on the source function, here we put 
the transverse flow velocity U r = by hand. In this case, the source function St{x,v) is 
almost proportional to the space-time volume of freeze-out hypersurface projected onto the 
transverse plane, and the azimuthal symmetry of the source function is restored as shown 
in Fig. [5[ As a result, the measure of the source opacity, (Ax) 2 — (Ay) 2 , vanishes. Although 
particle emission from the inside is less than that from the surface; the source is not opaque 
in this sense. When the transverse flow exists, the source function is deformed by the 
thermal Boltzmann factor exp(cosh.Yr cos (p Kt/Tj) and in our simulation the transverse 
flow rapidity, Yt, is almost proportional to r as shown in Fig. [2(]. Consequently, the number 
of emitted particles in x-direction is increased for the positive x region (cos0 > 0), especially 
at the surface, because of the behavior of Y T , and the number is decreased for the negative 
x region (cos0 < 0). This flow effect deforms surface dominant distribution in Fig. |5| to the 
crescent shape. This is the main reason why the source is opaque in our model. Though in 
the recent experimental result, i?o is larger than expected for the opaque source. According 



to the scenario as discussed in Ref. ||27|1 , the phase transition may cause a large time duration 
for the opaque source. However, in our hydrodynamical model, we assumed the first order 
phase transition, but the time duration At does not become large enough to reproduce the 
experiment. 

In order to discuss the effects of the phase transition on the source opacity, we sliced 
the source function St(x, y) by the freeze-out proper time r and show the source functions 
for each proper time in Fig. [|,[7|. Figures @(a)-0(j) show the same quantity as Fig. [| but 
they are not normalized. Integral N = / dxdy SV(x,y) should be considered as the number 
of particles. Figures |^(a)-|6|(f ) clearly show most particles are emitted from a thin surface. 
The measure of the source opacity (Ax) 2 — (Ay) 2 evaluated for Figs. H(b)-[7|(h) indicates 
extremely opaque property. For these sources, the YKP parametrization cannot be used 



because the YKP velocity becomes imaginary (i.e., A + B < 1C in Eq. (|22a|) ) [ft9 |. Figure 
^|(a) shows strong surface emission, but this source does not show the source opacity because 
the transverse flow is very weak at the early stage of evolution. (Ax) 2 — (Ay) 2 is negative 
throughout Fig. ||,[?|. As discussed in the Appendix, if first order phase transition exists, 
a flat-top structure appears in temperature distribution (Fig. 0) which corresponds to a 
large region caused by latent heat released at phase transition. This equi-temperature region 

9 



reaches freeze-out temperature almost simultaneously and makes a filled source function Fig. 
0(j). Hence, the existence of phase transition weakens the source opacity. 

V. RAPIDITY DEPENDENCE OF SOURCE PARAMETERS 



In this section, we discuss the pair rapidity Y nn dependence of YKP radii. The upper 
figure of Fig. |8| shows the R±, which is consistent with the experimental result. It is similar 
to the case of Kt dependence; our result of Rq shown in the lowest figure becomes smaller 
as Y n7T increases and our R does not agree with At. This departure also can be explained as 
the result of source opacity. The middle figure stands for R» , which shows a similar tendency 
as the experimental results. Source size Az is almost the same as the space-time extension 
Rn (|22d|) . Hence the interpretation in the YKP parametrization is applicable to our model. 
The disagreement of R» with the experiment is due to smaller longitudinal source width at 

larger Y n7T . In Fig. |§ we can see that the width of the source function Ar] = J (rf) — (r]} 2 
for //-direction becomes smaller as Y n7T increases. As shown in Fig.|TT|, this property does not 
appear in a Gaussian source model M where the source function is given as 



S.WK") 



rriT cosh(7/ — Y n ) 
(2tt) 3 ^2ti(6t) 2 



exp 



where 



G(rj) = exp 



T 



2{5r]) 2 



T 



TO) 



2(Sr) 



2R 2 



G(v), 



(25) 



(26) 



Of course, a small At] does not always mean a small Az in general, because z is given as 
z = rsinh r\. However, in the case where freeze-out takes place almost at some fixed r, Az 
is regarded as proportional to Ar\. And this case can be applied to our model, because for 
most of the particles, the freeze-out occurs at the latest stage (r ~ 10 or 11) as shown in Fig. 
|10| . The small width of the source function in our model is due to a non-Gaussian feature of 
^-dependence of a geometric factor of the source function, which corresponds to G(rj) in the 



Gaussian source model. Figure [12] shows the corresponding normalized geometric factor for 
the freeze-out hypersurface in our model; the geometric factor has a central plateau region. 
Non-Gaussian type source function comes from both initial temperature distribution and 
longitudinal flow. A deformed Boltzmann factor makes the source function narrower in a 
high Y n7T region. This feature is not only in a hydrodynamical model but also appears in 
thermal source models with a non-Gaussian geometric factor. In order to see the effect 
clearly, we replace G(t]) in Eq.(PB]) with the modified geometric factor G(r]). The solid line 
in Fig.|12| denotes the modified geometric factor G(r]) obtained by fitting the geometric factor 
for our model through 



G(rj) = exp 



(M -Vo)' 
2{5 V ) 2 



0(\v\-Vo) 



(27) 



Using the above G(rf) in Eq.(p5[) instead of G(t]), the same source function as in Fig. |IT 
is evaluated, and the result is shown in Fig. O. We can see that the width of the source 
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function becomes narrower at a high Y n7T region than a low Y nn region in the figure. Hence, 
we can say that the strong Y n7r dependence of R\\ comes from the geometric structure of the 
pion source in the rapidity space. 

Figure [14] shows pair rapidity Y nn dependence of the YKP rapidity and the longitudinal 
expansion velocity. In our calculation, we used Eq.(pO|) where Yykp is derived by using the 
YKP velocity v, which is given as 

F YK P = ^ln^ + F cm . (28) 

2 1 — v 

The dashed line stands for the infinite boost-invariant source. The departure between the 
solid line and the dotted line indicates that the failure of the simple interpretation of YKP 
velocity, which might be caused opacity of the source. The experimental results seem to 
be consistent with the infinite boost-invariant source except for the high rapidity region. 
Our results are much smaller than the boost-invariant model, as if our model contradicts 
the boost-invariance. Figure [15] shows average longitudinal rapidity (Yj,) as a function of 
longitudinal coordinate r\ for slow (Y nn = 3.4) particles and fast (Y n7T = 4.9) particles where 
Kt is fixed at 150 MeV. We can clearly see (Yjj) = r] for both cases. Figure [19] shows that 
Yl — i] is positive for small r but is negative in a later stage of the space-time evolution. 
Consequently, average rapidity (Yj,) is almost the same as 77. Hence there must be another 
reason why Yykp shows non-boost-invariant behavior. By making use of the source model, 
the result of the boost-invariant source is obtained by 5rj — > oo in Eg. ([2"5]) . For the finite 
but large enough Sr), the peak of the source function is only slightly shifted from r\ = Y nn , 
which is the peak of the Boltzmann factor exp(— my cosh (77 — Y n7T )/T). But for the small 
8r], a large shift will occur. Hence, the small Yykp reflects the finite size effect rather than 
longitudinal expansion. 



VI. CONCLUDING REMARKS 

In this paper, with the use of the numerical simulation for the CERN-SPS 158 A GeV 
Pb+Pb collisions based on a relativistic hydrodynamical model with first order phase transi- 
tion, we investigated the space-time structures via two-pion correlations in detail. In addition 
to the Cartesian parametrization, we adopted the Yano-Koonin-Podgoretskh parametriza- 
tion for two-particle correlation function and compared Kt and Y^ dependence of source 
parameters with the experimental result. In order to analyze the validity of the usual inter- 
pretations of HBT radii, we also evaluated space-time extent and source size for freeze-out 
hypersurface. 

In the case of Cartesian parametrization, we have compared Kt dependence of HBT 
radii and obtained results mostly consistent with the experiment. But in the case of YKP 
parametrization, Kt and Y n7T dependence of the transverse HBT radii, R±, are consistent 
with the experiment. The Kt dependence suggests the existence of the transverse flow, and 
the difference between our calculation and the experiment can be improved by including 
resonance decay into our model. The K T dependence of longitudinal HBT radii, R\\ , agrees 
with the experiment but Y^ dependence does not, which reflects a non-Gaussian structure 
of the geometric factor of the source function in //-space. For these two HBT radii, usual 
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interpretations of the transverse and longitudinal source size are available. But in the case 
of temporal source radii Rq, the interpretation of the time duration does not seem applicable 
at high Kt and Y nn . This comes from the fact that though first order phase transition exists, 
freeze-out hyper surface in our simulation is an opaque source. In the present stage, the 
experimental result is not clear whether the source is opaque or not |24] , |30"|| , and we need 
more analysis based on realistic source models. 
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APPENDIX: SPACE-TIME EVOLUTION 

Assuming that local equilibrium is achieved at 1 fm later than the collision instance, we 
put initial conditions of the hydrodynamical model on the r = r = 1 f m hypersurface. As 
for the initial local velocity, we use Bjorken's scaling solution J31|, Y\, = rj, in the longitudinal 



direction and neglect transverse flow, Y^ = 0. Initial entropy density distribution is given 



as, 



S(T ,rj) = S(T )exp 



iM ~" o) \M-m)-^fe(r-r ) 



£i ' \J y\ £i ' \J y 



(Al) 



and parameters we used in this paper are summarized in Table | . Those values are so 
designed to reproduce hadronic distributions of the recent experimental results with freeze- 



out temperature as Tf = 140 MeV. Figures 16 and 17 show the rapidity distribution and the 



transverse mass distribution of negative charged hadrons in Pb+Pb 158 A GeV collisions 
byNA49[0. 

Figure [18] displays the space time evolution of the temperature distribution of our hy- 
drodynamical model with the above parameters. From this simulation we can estimate 
"lifetimes" of fluid, QGP phase and mixed phase which are summarized in Table ||. 

In our simulation, the isothermal region at phase transition temperature T = T c appears 

very clearly. Though the global structure of the chronological evolution does not differ 

essentially from the model with smooth phase transition in a small temperature region 

15|,33], by virtue of the use of the volume fraction A at phase transition temperature, we 



can easily pick up the volume element with just the phase transition temperature. 

We describe the space time evolution of the local velocity in terms of Yi — rj, Yj>. Figure 



19| shows that Y^ — r\ becomes negative after r = 8.5 fm and the difference between our 



numerical solution and Bjorken's scaling solution is not significant on the whole. Figure [20 
indicates that the dependence of Yr on 77 is small and Yp is proportional to r. 

The temperature profile function (Fig. ^1]) shows the significant contribution of the space- 



time volume at T = T c . 
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TABLES 
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i'ABLE I. Initial Parameters 






T (MeV) 


Vo °r, r (fm) 




<7 r (fm) 


Pb + Pb 190 


0.7 0.7 1.2 x (207)( 1 / 3 ) - 


-1.0 


1.0 


TABLE II. Outputs 


Maximum initial energy density 






3.1 GeV/fm 3 


Maximum lifetime of fluid 






11.0 fm 


Maximum lifetime of QGP phase 






2.0 fm 


Maximum lifetime of mixed phase 






8.5 fm 
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FIG. 1. Kt dependence of Cartesian HBT radii at the rapidity of 3 < Y nn < 4. The reference 
frame is at Yq = 3.5. Open circles show the experimental results taken from NA49 preliminary 



data [ 23 1 . Closed circles show our results extracted from Eq. (18) through the Gaussian fitting. 
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FIG. 2. Kt dependence of YKP radii. Being the same as in Fig. ED, open and closed circles 
show the experimental data and HBT radii, respectively. Solid lines stand for space-time extensions 



|23a| )-(23c). Dotted lines stand for source sizes. (Az for Ru and At for Rq.) 
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FIG. 3. Source function St(x,v)- Thick regions emit particles more than thin regions. 
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FIG. 4. The second term -2(xi)/Px (dotted line) and third term ({Ax) 2 - (Ay) 2 )/Pj_ (solid 
line) in the right-hand side of Eq. (p3b|) are functions of Kt- Y n7r is fixed at 4.15. 
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FIG. 5. Same as Fig.H, but now the effect of transverse flow is neglected. 
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FIG. 6. Source functions St(x, y) sliced by freeze-out proper time. The thickness in each figure 
represents the intensity of the source. 
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FIG. 7. Rest of Fig. | 
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FIG. 8. Pair rapidity Y W7T dependence of YKP radii. As in Fig. y, open circles, closed circles, 
solid lines and dotted lines stand for experimental results, our HBT radii, space-time extensions 
and source sizes, respectively. 
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FIG. 9. Contour plot of source function for our model as a function of r\ and Y n7r . The window 
of Ynn is adjusted to Fig. @. 
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FIG. 10. The volume fraction of freeze-out hypersurface as a function of r. 
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FIG. 11. Contour plot of source function for the Gaussian source model (|25| ) as a function of 7/ 
and Ynn. 
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FIG. 12. The geometric factor of source function in 77-space. The box profiles stand for our 
hydrodynamical model. The solid line stands for modified geometric factor G{rf) obtained from 
the fit in Eq.@ where 770 = 0.92 and 5r] = 0.47 . 
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FIG. 13. Same as Fig. 11, but G{rj) is now replaced by G{rj) in the calculation. 
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FIG. 14. 1^-7,- dependence of Yykp- Open circles stand for the experimental results. Closed 
circles stand for our result obtained through Eq. ( |2§| ) . The solid line stands for our result calculated 
from space-time extensions Eq.( |22a| ). The dotted line stands for average longitudinal expansion 
rapidity (Yl) where (• • •) is defined by Eq.(p^). The dashed line stands for the infinite boost 
invariant source. 
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FIG. 15. Average longitudinal rapidity versus longitudinal coordinate r\. The thick line stands 
for Yttr = 4.9, and the thin line stands for Y^ = 3.4. The dotted line shows Bjorken's scaling 
solution Yl = rj. 
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FIG. 16. Rapidity distribution of primary negative charged hadrons (ir~, K~,p). Closed circles 




are taken from NA49 data 32]. The solid line stand for our numerical result. 
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FIG. 17. Transverse mass distribution of negative charged hadrons (n~ , K~ ,p) successively 
scaled down by an order of magnitude in the rapidity range y n £ {3.4,5.4}. Closed circles are 
taken from NA49 data [32]. The solid lines stands for our numerical results. 
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FIG. 18. The space-time evolution of temperature-distribution with initial temperature (190 
[MeV]), critical temperature (160 [MeV]) and freeze-out temperature (140 [MeV]). The flat area 
corresponds to the mixed phase region. From these graphs we can see that the QGP phase vanishes 
at r = 2.0 [fm] and the mixed phase vanishes at r = 8.5 [fm], and the temperature becomes less 
than the freeze-out temperature everywhere in the fluid at r = 11.0 [fm]. 
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FIG. 19. The space-time evolution of Yl — rj. 
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FIG. 20. The space-time evolution of Yt- 
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FIG. 21. The temperature profile function. 
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